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ABSTRACT 

For the origin of heavy r-process elements, different sources have been proposed, e.g., 
core-collapse supernovae or neutron star mergers. Old metal-poor stars carry the sig¬ 
nature of the astrophysical source(s). Among the elements dominantly made by the 
r-process, europium (Eu) is relatively easy to observe. In this work we simulate the 
evolution of europium in our galaxy with the inhomogeneous chemical evolution model 
’ICE’, and compare our results with spectroscopic observations. We test the most im¬ 
portant parameters affecting the chemical evolution of Eu: (a) for neutron star mergers 
the coalescence time scale of the merger (tcoai) and the probability to experience a 
neutron star merger event after two supernova explosions occurred and formed a dou¬ 
ble neutron star system (Pnsm) and (b) for the sub-class of magneto-rotationally 
driven supernovae (”Jet-SNe”), their occurrence rate compared to standard super¬ 
novae (Pjet-SN)- We find that the observed [Eu/Ee] pattern in the galaxy can be re¬ 
produced by a combination of neutron star mergers and magneto-rotationally driven 
supernovae as r-process sources. While neutron star mergers alone seem to set in at 
too high metallicities, Jet-SNe provide a cure for this deficiency at low metallicities. 
Eurthermore, we confirm that local inhomogeneities can explain the observed large 
spread in the europium abundances at low metallicities. We also predict the evolution 
of [0/Ee] to test whether the spread in a-elements for inhomogeneous models agrees 
with observations and whether this provides constraints on supernova explosion mod¬ 
els and their nucleosynthesis. 

Key words: Galaxy: abundances - Galaxy: evolution - nuclear reactions, nucleosyn¬ 
thesis, abundances 


1 INTRODUCTION 

The rapid neutron capture process (r-process, e.g., 
iThielemann et alll201ll . and references therein) is respon¬ 
sible for the production of about half of the heavy ele¬ 
ment abundances beyond Fe in the solar system, and of the 
heaviest elements like Th and U. The remaining heavy el¬ 
ement abundances are mostly made by the slow ne utron 
capture process (s-process, e.g.. iKappeler et aDl201ll) . De¬ 
spite its relevance, the true astrophysical origin of the r- 
process is still under debate. Because of the larger uncer¬ 
tainties affecting the r-process nucleosynthesis predictions 
compared to the s-process in stars, the r-process isotopic 
contribution in the solar system has been originally iden¬ 
tified by using the residual method, i.e., by subtracting 
the s-process component from the solar isotopic distribu- 
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tion (e.g.. lArlandini et al.lfl9^ . iBisterzo et Tll2014l) . The 
r-process residual abundances have been shown to be con¬ 
sistent in first approximation with the abundance signatures 
in old r-process rich metal-poor stars (at l east for elements 
heavier than Ba, see iTravaglio et ^12004 for details), car¬ 
rying the signatur e of the r-process nucleosynthes is in the 
early galaxy (e.g.. I^eden. Cowan fc Galling [20081 ) . For in¬ 
stance, Eu receives only a marginal contribution from the 
s-process (the s-process explains only 6 per cent of the so- 
lar Eu, while the re maining amount has an r-process origin, 
IBisterzo et al.l2014l) , and therefore it is often used as a tracer 
of the r-process nucleosynthesis in stellar spectroscopic ob¬ 
servations. One possibility to test predictions from r-process 
nucleosynthesis is to include the r-process stellar yields in 
galactic chemical evolution (GCE) simulations, and to com¬ 
pare the theoretical results with spectroscopic observations 
at different metallicities. Eu is an ideal diagnostic for these 
studies. The purpose of this work is to illustrate the eu- 
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ropium evolution throughout the evolution of our galaxy. We 
consider here the contribution from two sites (and their fre¬ 
quency) to the production of heavy r-process elements: Neu¬ 
tron Star Merger (NSM) and ’’magneto-rotationally driven 
supernovae” (hereafter referred to as ”Jet-SNe”). We will 
show that the combination of both sites is able to reproduce 
the observed europium abundance distribution of the stars 
of our galaxy. 

Neutron star - black hole (NS-BH) mergers might have 
a non-negligible contribution to the r-process inventory in 
the galaxy. However, their relevance as astrophysical source 
for the r-process is controversial, since thi s event has not yet 
been observed le.g.. lBauswein et al IB) . These difficulties 
also result in an extreme diver gence of the predicted galac - 
tic rate of such an event ('e.g.. IPostnov fc Yungelsonll2014f) . 
However, it should be noticed that a contri bution from NS- 
BH mergers has been predicted as wel l le.g.. lKorobkin et al.l 
l2012l . [Mennekens fc VanbevererJl2014l ~l. We give an estimate 
of the possible effects caused by these events in Section [S] 

Chemical evolution of galaxies has made strong ad¬ 
vances since its early days. Initially all approaches made use 
of the instantaneous recycling approximation in the sense, 
that the ejecta of stellar end stages were immediately uti¬ 
lized without delay after the Initial star formation, assuming 
that the stellar evolution time scale is short in compari¬ 
son to galactic evolution. If, in addition, the instantaneous 
mixing approximation was applied, i.e., assuming that the 
ejecta were instantaneously mixed throughout the galaxy, 
the whole galaxy acts as a homogeneous box. Neglecting this 
can explain radial gradients. Further developments included 
infall of primordial matter into and outflow of enriched ma¬ 
terial out of the galaxy (for a review of t hese early inves¬ 
tigations see e.g.. lAudouze fc Tinslevlll976l l. When relaxing 


the instantaneous recycling approximation, i.e., taking into 
account that (explosive) stellar ejecta enter the interstel¬ 
lar Medium (ISM) delayed with respect to the birth of a 
star by the duration of its stellar evolution, detailed predic¬ 
tions for the evolution of element abundances can be made. 
Based on nucleosynthesis predictions for stellar deaths, a 
number of detailed analyses have be en performed, from 
light elements up to the Fe-group le.g..iTimm es et al.||l9^. 


Goswami_&JPrantzo^j2000|jjMatteucd|_ 200 ll. Gibsonet_^lJ 


200j^ K^a yashl et al.ll2006l . iPagell l2009l . iKobavashil |2012| . 

Matteuc*^ 2012l l. Such approaches have recently also been 
applied to understand the enrichment of heavy elements in 
the galaxy (including r-process contrib utions) as a function 
of time or metallicity [Fe/H] (see e.g., Ishimaru_^_Wanajo| 
I 999 I. Travaglio et al'illOod. De Ponder fc Vanbeveren 20031 . 


Matteuccill2012l . IVangioni et al.ll2015l l. 

However, if still the instantaneous mixing approxima¬ 
tion is applied, i.e., such ejecta are instantaneously mixed 
with the global ISM, no local inhomogeneities can be pro¬ 
duced. The latter would relate to the fact that only limited 
amounts of the ISM are polluted by / mixed with the ejecta 
of each event. This effect is of essential importance, espe¬ 
cially at low metallicities, where portions of the ISM are 
already polluted by stellar winds and supernovae, and oth¬ 
ers are not. In addition, different portions of the ISM are 
polluted by different types of events, leading to a scatter at 
the same metallicity, which can in fact be utilized as a con¬ 
straint for these different stellar ejecta. When, however, uti¬ 
lizing the instantaneous mixing approximation, this leads to 


a unique relation between galactic evolution time and metal¬ 
licity [Fe/H], i.e., any [Fe/H] can be related to a specihc time 
in the evolution of a galaxy (while inhomogeneous mixing 
could experience similar [Fe/H] values in different locations 
of the galaxy at different times). This is especially the case 
in the very early galactic evolution ([Fe/H]< —2.5), when 
locally only a few stars (out of a whole initial mass func¬ 
tion IMF) might have exploded and imprinted their stellar 
neighbourhood with their ejecta. Thus, the application of 
chemical evolution models which utilize the instantaneous 
mixing approximation is questionable for the early evolu¬ 
tion of galaxies. 

In addition, for each [Fe/H], due to the instantaneous 
mixing, only a mean value of [X/Fe] (X being the element of 
Interest to follow in chemical evolution) is obtained. Inho¬ 
mogeneous mixing, however, could produce larger ratios in 
strongly polluted areas and smaller values in still less pol¬ 
luted ones. This means that the scatter in [X/Fe] at low 
metallicities, which might also be a helpful asset in point¬ 
ing to the origin of element X, cannot be reproduced or uti¬ 
lized with an homogeneous treatment. In case of rare events, 
which - on the other hand - produce large amounts of el¬ 
ement X in each event, this would produce a large scatter, 
and - if observed - could be used as a very helpful con¬ 
straint to identify the production site. For these reasons, 
specially for the origin of r-process elements like Eu, we 
think that only inhomogeneous chemical evolution models 
should be utilized at low metallicities. The two type of rare 
events (i) Jet-SNe (maybe up to 1% of all core collapse su¬ 
pernovae (CCSNe)) and (ii) neutron star mergers, with a 
similar occurrence frequency of about 1% of all CCSNe are 
considered here, while regular and more frequent core col¬ 
lapse supernovae might at most contribute to the lighter r- 
process elements. The bina r y mer ger rates a re estimated by 
lvaji_den Heuvel fc Lorirn^ Jl996ll as well as lKaloeera et ahl 
( 2004| . The rate of Jet-SNe is related to the fact that about 
1% of neutron stars are found with ma gnetic helds of the 
order 10^^ Gauss (magn etars, see, e.g., Kouveliotou et aP 
1 19981 . iKramer et al.llioo^ ') . 

Earlier inhomogeneous chem(odynam)ical evolution 
models for r-process elements like Eu have been provided 
by Travaglio et al. (2001, where the r-process yields were 
assumed to come fro m CCSN), lArgast et al.l (l2004l l and 
iMatteucci et al.1 (l2014ll comparing neutron star mergers and 
core collapse supernovae, Cescutti & Chiappini (2014, com¬ 
paring NSM and Jet-SNe) Mennekens & Vanbeveren (with 
NSM and NS-BH mergers), and Shen et al. and van de 
Voort et al. (2015, only utilizing neutron star mergers). 
One of the main questions here is related to the problem 
of reproducing [Eu/Fe] at low(est) metallicities. Cescutti & 
Chia ppini (2014) have sh own that this is possible with Jet- 
SNe. lArgast^et^^ (l2004h concluded that neutron star merg- 
ers cannot reproduce obs ervat ions at [Fe/H]< —2.5, while 
Ivan de Voort et al.1 (I 2 OI 5 I I and IShen et al.ll2015l came to the 
opposite conclusion. 

The main difference between Jet-SNe and neutron star 
mergers is that in one case the immediate progenitors are 
massive stars and the first occurrence in chemical evolution 
is due to the death of massive stars. In the other case the 
progenitors are also massive stars, leading to two supernova 
explosions in a binary system, which - if not disrupted - 
causes a binary neutron star system and a merger with a 
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given delay time due to gravitational radiation losses. Thus, 
one needs to consider two aspects: (i) the two supernova 
explosions and the pollution of the ISM with their ejecta 
(for the case of NS-BH mergers see the discussion in sec¬ 
tion El), and (ii) the delay time of the merger event after the 
formation of the binary neutron star system. Especially as¬ 
pect (i) can only be treated adequately with inhomogeneous 
evolution models, and there an additional factor is of ma¬ 
jor importance: with how much matter the supernova ejecta 
mix before the neutron star mergers eject their products into 
the same environment. 

This paper is organized as follows. In section [21 we in¬ 
troduce the model used to compute the evolution of abun¬ 
dances. In section[3l we present the influence of the different 
r- and non-r-process sites on the evolution. Additionally, we 
provide an overview why an inhomogeneous treatment of the 
evolution is important. In section|4lwe discuss the impact of 
inhomogeneities, causing and permitting a scatter of [X/Fe] 
ratios at low metallicities. As a further test of the model, 
we discuss the fact why the large scatter of [r/Fe] observed 
at low metallicities is strongly reduced for a-elements, and 
show how this constraints core collapse supernovae nucle¬ 
osynthesis predictions, which are still not available in a self- 
consistent way. In section [SI our results are summarized and 
discussed. 


Mtot 

Total infall mass 

10® Mq 

T 

time scale of infall decline 

5 ■ 10®yrs 

tmax 

time of the highest infall rate 

2 ■ 10®yrs 

t final 

duration of the simulation 

13.6 ■ 10®yrs 


Table 1. Main infall parameters. See Argast et al. (2004) for 
details on the parameters. 

steps of IMy. Primordial matter is assumed to fall into the 
simulation volume, obeying the form 

M(t) = a • t*’• (1) 

which permits an initially rising and eventually exponen¬ 
tially declining infall rate. While r and the total galaxy evo¬ 
lution time t final are fixed initially, the parameters a and b 
can be determined alternatively from Mtot (the total infall 
mass integrated over time), defined by 

final , 

Mtot .= a-t^ -e , (2) 

Jo 

and the time of maximal infall tmax, given by 

tmax :~ b • T . (3) 

See lArgast et al.l ll2004l for an extended discussion of 
the infall model and table [T] for the applied parameters. 


2 THE MODEL 


Rec ent chemodynamical galactic evolution models , like 
e.g., Minchev et ^ (l2014l ) nvan de Voort et all (l2015l) . and 
iShen et al.l I 20ll )7 can model in a self-consistent way mas¬ 
sive mergers of galactic subsystems (causing effects like in¬ 
fall in simpler models), energy feedback from stellar explo¬ 
sions (causing effects like outflows), radial migrations in disk 
galaxies, mixing and diffusion of matter/ISM, and the ini¬ 
tiation of star formation dependent on local conditions, re¬ 
sulting from the effects discussed above. In our present in¬ 
vestigation we still utilize a more classical approach with 
a pa rametrized infa ll of primordial matter, and a Schmidt 
law (ISchmidtII 19591 ) for star formation. Therefore, we ne¬ 
glect large scale mixing effects, while we include the feedback 
from stellar explosions and the resulting mixing with the sur¬ 
rounding ISM, according to a Sedov-Taylor blast wave. In 
this way, the model permits to keep track of the local in¬ 
homogeneities due to different CCSN ejecta. This approach 
allows to grasp the main features of the impact of the first 
stars / stellar deaths on the evolution of the heavy element 
enrichment. This approximation omits other mixing effects, 
e.g., spiral arm mixing (on time scales of the order of 2 • 10® 
years). The main focus of this work is the investigation of 
the chemical evolution behaviour at low metallicities, where 
these effects should not have occured, yet, and are therefore 
left out in this first order approximation. 

We treat the galactic chemical evolution of europium 
(Eu), iron (Fe) and Q-elements (e.g., oxygen O), utilizing 
the established GCE code ’’Inhomogeneous C hemical Evo¬ 
lution” (ICE), created by [Argast et al.l (|2004| ). A detailed 
description of the model can be found therein. 

For the simulation, we set up a cube of (2A:pc)® within 
the galaxy which is cut in 40® smaller cubes representing a 
(50pc)® sub cube each. The evolution is followed with time- 


2.1 Treating stellar births and deaths 

The main calculation loop at each time step (IMy) can be 
described in the following way. 

(i) We scan all mass cells of the total volume and calculate 
the star formation rate per volume and time step (10® yrs) 
acco rding to a Schmidt law with a density powe r a = 1.5 
('see lSchmidtlll95^ . lKennicuttlll99^ . lLarsorilll99ll ). Dividing 
by the average stellar mass of a Salpeter IMF (power —2.35) 
provides the total number of stars per time step n{t) created 
in the overall volume of our simulation. 

(ii) Individual cells in which stars are formed are selected 
randomly until n{t) is attained, but the probability is scaled 
with the density, which leads to the fact that patches of 
higher density, predominently close to supernova remnants, 
are chosen. 

(iii) The mass of a newly created star is chosen randomly 
in the range 0.1 to 50Mq, subject to the condition that 
the mass distribution of all stars follows a Salpeter IMF. 
Consequently only cells which contain more than 50Mq are 
selected in order to prevent a bias. 

(iv) The newly born star inherits the composition of the 
ISM out of which it is formed. 

(v) The age of each star is monitored, in order to deter¬ 
mine the end of its lifetime, either to form a white dwarf 
or experience a supernova explosion (see I2.2.TI and 1^.2.211 . 
A fraction of all high mass stars (M > 8Mq), according to 
the probability (Pjet-sv), is chosen to undergo a magneto- 
rotationally driven supernova event (see section 151^ . Type 
la supernova events are chosen from white dwarfs accord¬ 
ing to the discussion in l2.2.3l The treatment of neutron star 
mergers follows the description in 12.2.41 

(vi) The composition for the ejecta of all these events 
is chosen according to the discussion in 12.21 They will 
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pollute the neighbouring ISM with their nucleosynthe¬ 
sis products and sweep up the material in a chemically 
well mixed shell. We assume that an event pollutes typ¬ 
ically 5 • IO^Mq of surroun ding ISM due to a Sedov- 
Taylor blast-wave of 10^^ erg (iRvan, Norris fc Beers! Il996l . 
IShieevama fc Tsuiimotall998l ~l. This implies that the radius 
of a remnant depends strongly on the local density and the 
density of the surrounding cells. 

(vii) In the affected surrounding cells, stars are polluted 
by the matter of the previously exploded star and the event 
specific element yields. 

The details on the above procedure will be explained in 
the following. 


2.2 Nucleosynthesis sites 

2.2.1 Low (LMS) and intermediate mass stars (IMS) 

Low and intermediate mass stars provide a fundamental con¬ 
tribution to the GCE of e.g., He, C, N, F, Na and heavy s- 
process elements during the asymptotic giant branch (AGB) 
phase. For instance, most of the C a nd N in the solar sys¬ 
tem were made by AGB stars (e.g., iKobavashU I 2 OI 2 I 1 . In 
their hydrostatic burning phase, these stars lock-up a part 
of the overall mass and return most of it to the ISM in their 
AGB phase by stellar winds. Since the maximum radius of 
these winds is orders of magnitude smaller than the output 
range o f supernova events (e.g., radius of Grab remnant: 
5.5 Ly (|Hesteij[2008h . w hile the diameter of the Cat’s Eye 
Nebula is only 0.2 Ly (|Reed et al.lll999| H. our simulation 
assumes that stellar winds influence the ISM only in the 
local calculation cell. AGB stars provide only a marginal 
s-pro cess contribution to t ypical r-process elements like Eu 
fe.g.. pRavaglio et al.lll999h . In particular, for this work the 
s-process contribution to Eu plays a negligible role and we 
are not considering it here. 


2.2.2 High mass stars (HMS) 


Massive stars which exceed 8Mq are considered to 
end their life in a core-collapse supernova (CCSN, e.g.. 


Thielemann et al.lll996l. iNomoto et al.lll997l. IWooslev et al.l 

2 OO 2 I. Nomoto et al.l 20131. Jones et al.l 201311. CGSNe nro- 

duce most of the 
galaxy. They pre 
elements (S, Ca 
iron-group elem 
Sr neutron-mag 
ciated to GCSN 
components mif 
(e.s.. lArcones & 
possibly includi 
ular CGSNe as 
ments, as recen 
early hopes for 
right properties 
did not survive 
Liebendorfer et 

ments in the ini 

0 and Mg in the chemical inventory of the 
Dvide an important contribution to other a- 
Ti), to all intermediate-mass elements, the 
ents and to the s-process species up to the 
ic neak fe.g.. iRauscher et al.ll2002ll. Asso- 
e, different neutrino-driven nucleosynthesis 
fht be ejected and contribute to the GCE 
Thielemannll2013l. and references thereinl. 
ng the r-process. We did not include reg- 
a major source of heavy r-process ele- 
t investigations indicate strongly that the 
a high entropy neutrino wind with the 

(IWooslev et al. 

ll994l.lTakahashi et al.lll994ll 

advanced cor 
al.ll2003ll whicl 
lermost ejecta 

3 collapse simulations (e.g., 
led to proton-rich environ- 
[see alsolFischer et al.ll2010l. 
rather a so-called i/p-process 

HiideDohl et al.l 

2010lj. causine: 

(Frohlich et al. 

2006al, iFrohlich et alJ l2006bl. IPruet et alJ 


l2005l . IPruet et al.l l2006l . IWanaiol l2006h . Further investiga- 
tio ns seem to un derline this conclusion (recently revisited 
by IWanaidl2013l L although a more advanced - in medium 
- treatment of neutrons and protons in high density mat¬ 
ter cause s possible changes of the e l ectron fraction (Ye ) 
of ejecta (iMartinez-Pinedo et al.ir2012l : [Eberts et al.ll20lj ) 
and might permit a weak r-process, including small fractions 
of Eu. Similar effects m ight be possible via neutrino oscil¬ 
lations (IWu et al.ll2014l '). For this reason we did not include 
regular GCSNe in our GGE simulatio ns, although a weak r- 
process with small dHonda et ahllioo^ likel Eu contributions 
could be responsible for a lower bound of [Eu/Fe] observa¬ 
tions (see Fig. 0), explaining a non-detection of the lowest 
predicted [Eu/ Fe] ratios. Nucleos y nthes is y ields for HMS 
are ta ken from [Thielemann et ^ (Il996l l or iNomoto et al.l 
lll997tl . Assuming a typical explosion energy of 10®^erg, the 
ejecta are mixed with the surrounding interstellar medium 
via the expansion of a Sedov-Taylor blast wave, which stops 
at a radius which contains about 5 • IO'^Mq (see section 12.11 
for details on the iteration procedure). 


2.2.3 Supernovae Type la (SNIa) 


When an IMS is newly born in a binary system, there is 
a probability that it has a companion in the appropriate 
mass range leading finally to a SNIa, following a double- or 
single degenerat e scena rio. We follow the analytical sugges¬ 
tion of lGreggi3 (l2005ll and reduce the numerous degener¬ 
acy parameters to one probability (PsNia = 9 • 10~^) for a 
newly born IMS to actually be born in a system fulfilling 
the prerequisites for a SNIa. Once the star enters its red gi- 


ant phase, we let the sy 
and emit the event spe 
model CDD2), which h 
with iron. For this work 
metallicity, consistently 
lations. We are aware t 
several SNIa yields inch 
metallicity of the SNIa 

stem perform a SNIa-type explosion 
cific vields fcf. Ilwamoto et al.lll999l. 
ighly enriches the surrounding ISM 
we use the same SNIa yields for each 
with the|Argast et al.l (l2004|j calcu- 

hat this choice is not optimal, since 
rding e.g., Mn and Fe depend on the 
progenitor fe.g..lTimmes et al.ll2003l. 

Thielemann et al.l 2004 

. Trayaglio et al. 20051. iBrayo et al.l 

2 OIOI. Seitenzahl et al. l2013h. On the other hand, this an- 


proximation does not have any impact on our analysis and 
our conclusions. 


2 . 2.4 Neutron Star Merger (NSM) 


If two newly born HMS were created in a binary system, 
they may both undergo a GCSN individually. This could 
leave two gravitationally bound Neutron Stars (”NS”) be¬ 
hind. Such a system emits gravitational waves and the 
two NS spiral inwards towards their common center of 
mass with a coalescence time (tcoai) until they merge. 
The actual merging event is accompanied by an ejection 
of matter and (r-process) nucleosynthesis |Rossto 1 12 OI 3 I. 
FLeiburglmus_et_aL[|l999yPanov_e£^ I 200 R [Korobldr]_et_^ 


2012, Sauswein_et_alJ_l20W j_ Rosswog 2oS iRosswog et af. 

2014 , lEichler et al.ll2014l . lWanaio et al. 120141 . As all of these 


publications show the emergence of a strong r-process, in 
the mass region of Eu they suffer partially from nuclear 
unce rtainties related to fission fragment dis tributions (see 
e.g.. lEichler et al]|2014l . iGorielv et al.l l2013fl . For our pur¬ 
poses we chose to utilize as total amount of r-process ejecta 
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1.28 • 10 ~^M(r) (consisten t with the 1.4M(7) + 1.4M0 NS col¬ 
lision in iKorobkin et al.lI 2 OI 2 I and iRosswoj I 2 OI 3 I I, but dis¬ 
tributed in solar r-process proportions, which leads for Eu 
to a total amount of 10 “^Mq per merger. This value is rel¬ 
atively high in comparison to other investigations in the lit¬ 
erature. 

Observational constraints for the probability of a newly 
born sta r to undergo th is proc edure (Pjvsm) are provided 
by e.g., iKalogera et al.l (|2004l ') who have found a NSM 
rate of Rnsm = 83.0lgg®j^Myr“^, which corresponds to a 
Pnsm = 0.0180jlQ Qj 43 . The coalescence time, Pnsm and 
the event specific yields are important parameters for GCE, 
and their influence on the GCE are subject of this pa¬ 
per. Concerning the coalescence time scale, it might be 
more realistic to use a distributio n function (e.g., as in 
llshimaru. Wanaio fc Prantzoj2015l i instead of a fixed value. 
We utilize this simplified procedure as a first order approach. 


2.2.5 Magnetorotationally driven supemovae (Jet-SNe) 

A fraction (Pjet-SN) of high mass stars end their life 
as a "magneto-rotationally driven supernova” or mag- 
netar, forming in the center a highly magnetized neu¬ 
tron star (with fields of the order 10^® Gauss) and eject- 
ing r-process matter along the poles of the rotation axis 
dFuiirnoto et al.ll200fil. iFuiimoto et ahll^OOSl : IWinteler et al.l 
I 2 OI 2 I . Mosta et al.l I 2 OI 4 I ). r-process simulations for suc h 
events were first undertaken in 3D bv lWinteler et ahl (|2ni2h . 
For the purpose of this work, we randomly choose newly 
born high mass stars to later form a Jet-SN. At the end 
of their life time, they explod e similar to a GCSN, however 
with different ejecta. Based on IWinteler et aP (|2012h . we as¬ 
sume an amount of 14 • 1O“®M0 of europium ejected to the 
ISM by such an event. In this work, we study the influence 
of Pjet-SN and the specific Jet-SN yields on the GCE. 


2.3 Observed stellar abundances 


Data for the observed stars to compare our simulation 
results with are taken from the SAGA (Stellar Abundances 
for Galactic Archaeo l ogy) database (e.g ., ISuda et al.1 
l2008l . ISuda et al.1 1201 ll . lYamada et al.l 1201 ^ in particular 


[Eu/Fe1 abundance observations are ma i nly from e.g. 


ftancois et aJj 200?!. [simmerer et al.l 


2004. Barklem et al 


2014al. Roederer et al 

|2014bl, 

Roederer et al.1 l2014d. 

Shetrone. Gote. Stetson 

2001. 
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From the raw data, we excluded carbon enriched metal 
poor stars (’’GEMPs”) and stars with binary nature, since 
the surface abundances of such objects are expected to 
be affected by internal pollution from deeper layers or 
pollution from the binary companion. 


3.1 Coalescence time scale and NSM probability 

As a prerequisite, we studied the influence of both coales¬ 
cence time and the probability of a binary system to become 
a NSM. In Figure [T] we present the evolution of [Eu/Fe] 
abundances when only NSM contribute to the enrichment. 
The results can be summarized as follows. 

(i) Smaller coalescence time scale leads to an enrichment 
of europium at lower metallicities. Larger coalescence time 
scale shifts this to higher metallicities. 

(ii) A higher NSM probability shifts towards a quantita¬ 
tively higher enrichment combined with an appearence at 
lower metallicities. 

These effects can be explained in the following way. 

(i) When binary neutron star systems take longer to co¬ 
alesce, the time between the GCSN of both stars and the 
NSM event is longer. The longer this delay time, also fur¬ 
ther nucleosynthesis events occur in the galaxy during this 
period, enriching the ISM with metals. Thus, when the NSM 
event finally takes place, surrounding stars have developed a 
higher [Fe/H] abundance, shifting the system towards higher 
[Fe/H] abundances, respectively. This implies an overall eu¬ 
ropium production shift towards higher metallicities. 

(ii) With more binary systems becoming NSM, the pro¬ 
duced europium amount per time step is larger, since every 
event produces the same amount of r-process elements. This 
leads to a higher [Eu/Fe] abundance, compared to simula¬ 
tions with lower NSM probability. As the fraction of NSM 
systems is higher while the GCSN rate is constant, larger 
amounts of europium are produced, while the surrounding 
medium evolves regularly. This also leads to a higher abun¬ 
dance of europium at lower [Fe/Hj. These effects shift the 
[Eu/Fe] curve to higher values for the same [Fe/H[. 

All these results are co nsistent with the earlier conclusions 
by lArgast et al.l (120041 '). stating that it is extremely diffi¬ 
cult to reproduce the observed [Eu/Fe] ratios at metallicities 
[Fe/H]< —2.5 by NSM alone. A potential solution would be 
that the preceding supernovae which produced the two neu¬ 
tron stars of the merging system mix their ejecta with more 
extended amounts of the ISM. We utilized the results fol¬ 
lowing a Sedov-Taylor blast wave of 10®^ erg, which pollutes 
of the order 5 ■ 1O" ^ M0 o f ISM until the shock is stopped. 
Ivan de Voort et al.l (I 2 OI 5 I I assumed (in the ir standard case) 
the mixing with more than 1O®M0 of ISM JShen et al.ll201a 
utilized 2 • 10® M 0 in a similar approach). This produces an 
environment with a substantially lower [Fe/H] into which 
the NSM ejecta enter. Thus, it is not surprising that in such 
a case the Eu enrichment by NSM is setting in at lower 
metallicities. The higher resolution run shown in Fig. 4 of 
Ivan de Voort et al. I (l20l5l) agrees with our results. Thus, the 
major question is whether such a very much enlarged mix¬ 
ing with the ISM by almost two orders of magnitude can be 
substantiated. We will discuss these aspect further in sec¬ 
tion [S] 


3 RESULTS 

For a general understanding of the effects of Jet-SNe and 
NSM on GGE, namely the parameters Pnsm, tcoai and 
PJet-SN, we performed a number of simulations described 
in detail below. 


3.2 Probability of Jet-SNe 

The contribution of Jet-SNe to the galactic Eu abundance 
differs from that of NSM. Since Jet-SNe explode directly 
from a massive star, they contribute much earlier to the 
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Figure 1. Influence of coalescence time scale and NSM- 
probability on Eu-Abundances in GCE. Magenta stars represent 
observations. Red dots correspond to model star abundances as 
in Argast et al. (2004). The coalescence time scale of this event 
is 10® years and the probability PjvSM is set to 4T0“^. Green 
dots illustrate the effect on the abundances if the coalescence 
time scale of NSM is shorter (around 10® years). Blue dots show 
the abundance change if the probability of HMS binaries to later 
merge in a NSM is increased to 4T0“^ (cf. subsection [TTJ . 


Figure 2. Influence of increased Jet-SN probabilities on Eu- 
Abundances in GCE. Magenta stars represent observations. 
Green dots represent model star abundances based on Winteler 
et al. (2012), the Jet-SN probability has been chosen to follow the 
observations at [Ee/H]> —1.5. A good value seems to be 0.1% of 
HMS to end up in a Jet-SN. Note that this model fails to repro¬ 
duce the observed abundances at lower metallicities. Blue dots 
illustrate the effect on the abundances if the Jet-SN probability 
is increased to 1%. This model better reproduces the observed 
abundances at lower metallicities, but clearly fails at higher ones. 


chemical evolution than NSM. Since the interstellar matter 
is distributed more inhomogeneously than in later evolution 
stages of the galaxy, high [Eu/Fe] abundances are possible 
in individual stars. This leads to a large spread in the abun¬ 
dances towards lower metallicities. Considering Jet-SNe, the 
parameter with the highest impact on GCE for such rare 
events, similar to NSM events but ’’earlier” in metallicity, is 
the probability of a massive star to actually become a Jet- 
SN. A lower probability leads to a smaller overall [Eu/Fe] 
abundance, while a higher probability leads to larger abun¬ 
dances. However, we also recognize a larger spread in abun¬ 
dances in models with lower probability. This comes from 
the fact that the high yield of the event only sets an upper 
limit on the abundances. The rarer an event is, the more and 
the longer stars remain unpolluted. This results in a larger 
spectrum of abundances in stars and therefore in a larger 
spread in [Eu/Fe] ratios. Note from Fig. [2] and Fig. [3] that 
Jet-SNe might explain the abundances at low metallicities 
better than NSM. Thus, while Jet-SNe alone could be an 
explanation for the lower metallicity observations, there is 
clear evidence of NSM events and therefore we have to exam¬ 
ine the combination of both events. Whether the apparently 
to high concentration of model stars with low [Eu/Fe] values 
at metallicities —3 <[Fe/H]< 2 in comparison to observa¬ 
tions is related to observational bias or whether we require 
another additional source will be discussed in the following 
sections. 

3.3 Combination of sites 

If both sites (Jet-SN and NSM) are considered to contribute 
to the galactic europium abundances, their contributions 
overlap. Therefore, parameters which lead to the observed 
[Eu/Fe] abundances, have to be searched for. As described 
in section [331 NSM contribute at a delayed stage to the 
GCE and in our simulations are unable to reproduce eu¬ 
ropium abundances at metallicities [Fe/H]< —2.5, Jet-SNe, 
however, contribute europium early, but only in those re- 
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Figure 3. Same as Figure [2] but with decreased probabilities. 
Red dots are the same as green dots in Fig. [2] with Jet-SN proba¬ 
bility of 0.1%; Green and blue dots represent a Jet-SN probability 
of 10“"^ and 2 • 10“^, respectively. From the comparison of these 
models, we can see how decreased event probability shifts the 
abundance curve down. We also remark an increase of the spread 
in abundances when the probability is lowered. The rarer a high 
yield event is, the larger is the spread in abundances. 

gions where they occured, and cause a larger spread in the 
[Eu/Fe] values (cf. section We have to test whether it 
is possible to use the same parameters as in sections [3^ and 
El since the full combination of both events could lead to an 
overproduction of elements. We can use the earlier param¬ 
eter explorations to tune the simulated abundance pattern 
in order to match the observations. In the following, we will 
discuss two possible cases: 

(i) Pnsm = 3.4 • 10“"*, Pjst-SN = 0.3%, tcoai = IMy 
(hereafter model Jet-|-NSM:A). The results for the model 
Jet-|-NSM:A in comparison with observations are shown in 
Figure (H This model provide a reasonable explanation of 
the observations at lower and higher metallicities, but there 
is an overproduction of europium between —2 < [Fe/H] < 
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Figure 4. Evolution of Eu-abundances in GCE including both 
Jet-SNe and NSM as r-process sites. Magenta stars represent 
observations, whereas blue dots represent model stars. Model 
(Jet+NSM:A) parameters are Pnsm = 3.4 ■ 10“^, Pjet—SN = 
0.3%, tcoai = IMy. 
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Figure 5. Same as Figure [J] but with a different parameter set 
(ModelJet+NSM:B). Magenta stars represent observations (with 
observational errors; however, magenta stars at low metallicities 
which carry only horizontal errors represent upper limits). Blue 
dots represent model stars with Pnsm = 3.8 • 10“"^, Pj^t—SN = 
0.1%, t,,al = lOMy. 

— 1. We conclude that larger coalescence time scales and 
larger probabilities are necessary regarding NSM, and lower 
probability of Jet-SNe is necessary to flatten and lower the 
modelled abundance curve. 

(ii) Pnsm = 3.8 • 10“"*, Pjet -SN = 0.1%, tcoal = lOMy 
(Model Jet-|-NSM:B). The results for the model Jet-|-NSM:B 
in comparison with observations are shown in Figure O This 
model explains the main features of the abundance curve 
quite well: The spread at low metallicities, the first confine¬ 
ment of the spread at [Fe/H]Ri — 2, the plateau between 
[Fe/H]« — 2 and [Fe/H]« — 0.6, and the second confinement 
of the spread at [Fe/H]« —0.2. However, there still seem 
to be difficulties at [Fe/HJps — 2: the scatter in abundances 
towards low [Fe/H] ratios seems to be a bit too broad. This 
spread might be slightly reduced by additional mixing terms 
(e.g., spiral arms mixing) or an additional source providing 
ratios of [Eu/Fe]= —1, which we did not consider in this 
work. 

Considering Fig. |4] and Fig. [5] While the results from 
both models Jet-|-NSM:A and Jet-|-NSM:B can reproduce 



[ Fe/H ] 


Figure 6. Enrichment history tor models Jet+NSM:A and 
Jet-|-NSM:B (cf. Figure |4] and [S] for evolution plots). Magenta 
stars represent observations, whereas blue dots represent model 
stars as per Figure[5] (Model Jet-|-NSM:B). Red dots representing 
the enrichment history of the simulation as per Figure (Model 
Jet-|-NSM:A) do not suit the observational data. 

the observed spread of [Eu/Fe] in the early galaxy, model 
Jet-|-NSB:B seems to better fit the overall [Eu/Fe] vs. [Fe/H] 
distribution. On the other hand, the evolution of the [Eu/Fe] 
ratio at low metallicity depends on the r-process produc¬ 
tion and on the Fe production in CCSNe (see Section [4] 
and discussion), In Fig. (6] we compare the results for the 
enrichment history of europium in the galaxy according 
to Jet+NSM:A and Jet-|-NSM:B models with observations. 
While the [Eu/H] vs. [Fe/H] ratios predicted by model 
Jet-|-NSM:B are in agreement with the observations, model 
Jet-|-NSM:A seems to be ruled out. 

4 THE IMPORTANCE OF 
INHOMOGENEITIES 

4.1 Inhomogeneities in GCE 

From observations of [Eu/Fe] in the early galaxy, one of 
the main features is a spread in the abundance ratios. Our 
model is able to reproduce these spreads, mainly because 
of the inhomogeneous pollution of matter. In Fig. [7] we try 
to illustrate the effect of applying such an inhomogeneous 
model. For this purpose, we perform a cut through the xy- 
plane of the simulation volume for specific time steps. These 
time steps are marked in the top panel of Fig. [3 in order 
to provide the reader with a quick glance of the extent of 
the inhomogeneous element distribution at the correspon¬ 
dent metallicities. For each marker, we provide the complete 
density field at this specific time step in the middle and 
lower panels (cf. figure caption for details). We show the 
extent of inhomogeneities in the middle left panel, for the 
first marker in the upper panel of the Figure. At this time 
step, we can see - by counting the ’’bubble”-style patterns - 
that at least three supernovae must have taken place before 
the snapshot. Since such events give rise to nucleosynthesis, 
the abundances of metals in such a supernova remnant bub¬ 
ble are higher than outside such a remnant. A star being 
born inside such a remnant will inherit more metals than 
a star born outside. Therefore, in the early stages of galac¬ 
tic evolution the stellar abundances are strongly affected by 
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the location where a star is born. Considering much later 
stages of the evolution, (e.g., the lower right panel of Fig. 0 
corresponding to the fourth marker of the upper panel) the 
supernova remnants have a large overlap. Numerous super¬ 
nova explosions, have contributed lots of nucleosynthesis all 
over the galaxy. This leads to an averaged distribution of 
abundances, including different events and an integral over 
the initial mass function of stars. Therefore, it resembles 
a ’’mixed” phase of galactic evolution, where the elements 
have been homogenised over the whole volume. At this stage 
of the evolution, it seems not to be so relevant where a star 
was born. As a consequence, there are smaller differences 
in the abundance of metals in stars. Therefore, a confine¬ 
ment in the spread of abundances of chemical elements at 
later stages of the chemical evolution is obtained. Becoming 
more and more homogeneous, the [Eu/Fe] value converges to 
a value that can be obtained by integrating the event yields 
over the whole IMF. 

4.2 Instantaneous Mixing Approximation 

A number of recent chemical evolution studies revoked 
the ’’instantane o us m ix ing approximation” (I.M.A., e.g., 
Chfappini et al.l 1200 ll . iRecchi et al.l l200ll . ISpitoni et alJ 
2 OO 9 II . The I.M.A. simplifies a chemical evolution model in 
terms of mass movement. In detail, all event outputs are 
expected to mix with the surrounding ISM instantaneously. 
Such approaches always result in an average value of ele¬ 
ment ratios for each [Fe/H]. Therefore, the I.M.A. scheme 
all stars at a given time inherit the same abundance patterns 
of elements and it is impossible to reproduce a scatter in the 
galactic abundances, which seems to be a crucial ingredient 
at low metalliclties. Indeed, instead of a spread of distri¬ 
butions only one value is obtained for each metallicity. We 
calculate the best £t model (Jet-|-NSM:B, cf. Fig. [5]) with 
I.M.A. The result can be found in Fig. [8] The I.M.A. ap¬ 
proach may be used to study the chemical evolution trends 
with a lower computational effort, but Figure [3 shows that 
the reproduction of spreads in abundance ratios due to local 
inhomogeneities requires to use more complex codes as e.g., 
the ICE code adopted for this work. 

While inhomogeneus GCE codes can explain the spread 
in r-process elements, there is the question whether they 
might predict a far too large spread for other elements 
(e.g., a elements) at low metalliclties (with present stel¬ 
lar yields from artihcially induced CCS N explosion models). 
Such e ffects can also be seen in Fig. 1 in Ivan de Voort et al.l 
ll2015h for [Mg/Fe], spreading by more than 1 dex, while ob¬ 
servations seem to show a smaller spread up to 0.5 dex. 
This can be related to the amount of supernova ejecta 
being mixed with the ISM (see discussion above and in 
section [S] a more extended mixing reduces this spread), 
but it can also be related to the supernova nucleosynthesis 
yields which were never tested before in such inhomogeneous 
GCE studies. From general considerations of chemical evo¬ 
lution studies, it is found that there are large uncertainties 
for G CE studies, pa r ticula rly the influence of stellar yields 
(e.g., iRomano et al. I (|20l^ ). In Fig. H we show the results 
of model Jet-|-NSM:B, u sing the CCSN yield s from from 
iNomoto et al.l lll997ll and iNomoto et ^ l|2006 ll. which con- 
hrms a large spread in [0/Fe], similar to Ivan de Voort et al.l 
ll2015h for [Mg/Fe]. However, present supernova yields are 


Figure 7. The top panel shows the same GCE-model as in Fig. [5] 
(Model Jet-|-NSM:B), but without observations; The red markers 
refer to the position where a density determination cut through 
the the xy-plane of the simulation volume is performed. The mid¬ 
dle and lower two panels show the density distribution through 
these planes. The middle left panel corresponds to the very left 
marker (”A”) position’s density profile (approximately 180 mil¬ 
lion years (My) have passed in the simulation), the middle right 
panel to the second marker ”B” (r^ 290 My), the lower left panel 
to the third marker ”C” (~ 2 Gy) and the lower right panel to 
the very right marker ”D” (si 12 Gy). 

the result of artificially induced explosions with constant 
explosion energies of the order of 10®^ erg. If we consider 
that explosion energies might increase with the compactness 
of th e stellar core (i.e., progenitor mass, e.g., IPereso et al.l 
l2015h . the heavier a-elements and Fe might be enhanced as 
a function of progenitor mass. On the other hand O, Ne, 
and Mg yields are dominated by hydros tatic burning and 
also increase with progenitor mass (e.g., iThielemann et aH 
Il996l) . This could permit to obtain more constant a/Fe ra¬ 
tios over a wide mass range, although the total amount of 
ejecta differs (increases) as a function of progenitor mass. 
This scenario does not take into account all the complexity 
and the multi-d imensional nature of the CCSN event (e.g., 
iHix et ^ I I 2 OI 4 and references therein) that should be con¬ 
sidered, but it may be interesting to test its impact in our 
GCE simulations. In Figs. [9land llOI we show the results for 
tests where we: 
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Figure 8. Same GCE-model as in Fig. [5] (Model Jet+NSM:B); 
however I.M.A. is applied instead of inhomogeneous evolution. 
One is able to observe a trend in the abundance evolution, how¬ 
ever the scatter in the abundance pattern is not present any¬ 
more (cf. Section |4] for further discussion). The kink at about 
[Fe/H]= —2.5 is related to the delayed time when NSMs set in 
and contribute to Fu as well. This Figure can also be compared 
to Fig. 2 in Matteucci et ah (2014) which shows the contribution 
of NSMs alone for various merger delay times and Eu production 
yields and Fig. 5 in Vangioni et ah (2015) [mergers alone being 
indicated by black lines]. Thus, also in this approach it is evident 
that the explanation of [Eu/Fe] at low(est) metallicities by NSM 
alone is not possible. 


(i) replace the lNomoto et al] (1 19971 ) iron yields by ad-hoc 
yields, fitting, however, the observed SN1987A iron produc¬ 
tion; 

(ii) keep the same CCSN rate as in the previous models; 

(hi) adopt the parameters to study the r-process nucle¬ 
osynthesis of Model Jet-FNSM:B, obtaining the same [Eu/H] 
ratio. 

This leads, based on the adopted CCSN yields, to a possi¬ 
bility to minimize the spread in a-elements at low metal¬ 
licities, while keeping the spread in the r-process element 
evolution. Therefore, the spread of [O/Fe] obtained from 
GCE simulations for the early galaxy is strongly affected 
by the uncertainties in the stellar yields, and it is difhcult 
to disentangle them from more intrinsic GCE uncertainties. 
This means that at this stage it is not obvious whether an 
overestimation of the observed [O/Fe] spread is a problem of 
the ICE code, the observations could rather provide a con¬ 
straint on stellar yields. In particular, the use of realistic, 
self-consistent, explosion energies, might reduce the spread 
at low metallicities to a large extent. Another fundamental 
point is related to the discussion in Section [3] concerning 
[Eu/Fe]. At this stage, we consider [Eu/H] as more con¬ 
straining to study the r-process nucleosynthesis compared to 
[Eu/Fe], since Fe yields from CCSNe are affected by large 
uncertainties. Therefore, the model Jet-|-NSM:B is recom¬ 
mended compared to Jet-|-NSM:A (see also Fig. 1^. 


5 CONCLUSION AND DISCUSSION 

The main goal of this paper was to reproduce the solar eu¬ 
ropium abundance as well as the evolution of [Eu/Fe] vs. 
[Fe/H] throughout the evolution of the galaxy. For this rea- 


Figure 9. Same GCE-model as in Fig. [5] (Model Jet-|-NSM:B); 
Red dots show the abundance evolution of Oxygen when Nomoto 
et al. (1997) yields are employed, while blue dots represent a far 
narrower spread at low metallicities if ad hoc yields are applied 
(which still would need to be optimized to obtain a better agree¬ 
ment with the metallicity evolution between — 1 < [Fe/H] < 0). 
Note that the downturn at high metallicities is shifted to higher 
[Fe/H] values. This is probably due to an overestimate of the total 
IMF-integrated Fe-production, which should be improved with re¬ 
alistic self-consistent explosion models and their iron yields. While 
the delay time scale for SNIa is unchanged, earlier CCSN pro¬ 
duce more iron, thus dispersing the whole abundance curve. Here 
we only want to show how changes to possibly more realistic, 
progenitor-mass dependent, explosion energies can improve the 
[a/Fe] spread, while the [r/Fe] spread is conserved. (Cf. Section[4] 
for further discussion.) 
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Figure 10. Same consideration as in Fig. [9] however with 
magnesium instead of oxygen. GCE-model as in Fig. [5] (Model 
Jet-|-NSM:B); Red dots show the abundance evolution of magne¬ 
sium when Thielemann et al. (1996) / Nomoto et al. (1997) yields 
are employed, while blue dots represent a far narrower spread at 
low metallicities if ad hoc yields are applied. Cf. Section [4| for 
further discussion.) 



son we have studied the influence of two main r-process sites 
(NSM and Jet-SNe) on the GCE. 

Our simulations were base d on the inhomoge neous 
chemical evolution (ICE) model of lArgast et al.l (12003 ). with 
updated nucleosynthesis input for the two sites considered, 
their respective occurrence frequencies / time delays, and a 
model resolution of (50pc)^. The main conclusions are that: 

(i) The production of heavy r -process matter in NSM 
is evident since many years (see iFreiburghaus et al.l Il999l 
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Figure 11. Effect of slightly increased sweep-up mass on GCE. 
Magenta stars represent observations. Red dots show model stars 
as per our reference model JET-|-NSM:B. Blue dots represent a 
model where every CCSN pollutes 2 ■ 10 ®Mq of ISM. The dom¬ 
inant effect of this increased sweep up mass is to decrease the 
scatter of abundances and to shift the abundance curve towards 
lower metallicities.) 


and many later investigations ud to iKorobkin et al.l 

20121, 

R.osswod 2 OI 3 I. iBauswein et al.l l2013l. iR.osswoe et al.1 

2014, 

Just et al.l I 2 OI 4 I. IWanaio et al.l l2014l. Eichler et al.l 

2014, 

Mendoza et al.l 2015|). Our implementation of NSM in the in- 


homogeneous chemical evolution model ”ICE” can explain 
the bulk of Eu (r-process) contributions in the galaxy for 
[Fe/H]> —2.5, but have problems to explain the amount and 
the spread of [Eu/Fe] at lower metallicities. This is in agree¬ 
ment with the initial fi ndings of Argast et al.| | 20p4ll . Recent 
SPH-based studies by Ivan de Voort et al. ( 2015h make use 
of a mixing o f the e jecta with 3 • 1O®M0, a further study 
bv IShen et al.l (l2015l l utilizes a mixing with 2 • 10® Mq up to 
8T0 ®Mq. The mixing volume we utilize, based on the Sedov- 
Taylor blast wave approach, would be related to a subgrid- 
resolution in these studies, but this treatment is essential 
for the outcome. Mixing initially with a larger amount of 
matter causes smaller [Fe/H] ratios into which the r-process 
material is injected. 

We have tested such differences in mixing volumes/masses 
also within our ICE approach. Fig. If II shows the results we 
obtain when changing from the Sedov-Taylor bla st wave ap¬ 
proac h to a mixing mass of 2 ■ 10® Mq (like in IShen et al.l 
I 2 OI 5 I I. and we can see that we essentially reproduce their re¬ 
sults. On the other hand, a high er resolution test in section 
3.1 of ivan de Voort et al.l (l2015l ~l is essentially in agreement 
with our results presented above. Thus, these differences are 
not based on the differences in sophistication of the multi-D 
hydrodynamics approach, permitting to model energy feed¬ 
back from supernovae, outflows and infall, they can rather 
be linked directly to the mixing volumes of supernova ejecta. 
This requires further studies in order to understand whether 
there exist physical processes (on the timescale of the de¬ 
lay between the supernova explosions and the merger event) 
which permit a mixing beyond the Sedov-Taylor blast wave 
approach. 

(ii) The production of heavy r-process elements in a rare 
species of CCSNe with fast rotation rates and high mag¬ 
netic fields, causing (fast) jet ejection of ne utron-rich matte r 
along the poles has first been postulated bv ICamero nl ll200lli . 
This was followed up in rotationally symmetric 2D calcula¬ 


tions [^jimoto^^F^ai] ll2006l~l. F uiimoto et al.l (l 2008l 'l and the 
first 3D calculations by IWinteler et ~ I 2012ll . These cal¬ 
culations still depend on unknown rotation velocities, and 
magnetic field configurations before collapse, however, they 
agree with the observations of magnetars and neutron stars 
with magnetic fields of the order 10^® Gauss which make up 
about 1 % of all obser v ed ne utron stars. Further 3D calcula- 
tions by Mosta et al.| ( 2014ll and recent 2D calculations by 
iNishimura et akn 2015 1. might indicate that not all events 
leading to such highly magnetized neutron stars are able to 
eject the heaviest r-process elements in solar proportions. 
Thus, probably less than 1% of all CCSNe end as Jet-SNe 
with a full r-process. 

When introducing Jet-SNe with ejecta as predicted by 
IWinteler et al.l (I 2 OI 2 I I. they can fill in the missing Eu at 
lower metallicities and reproduce the spread in [Eu/Fe1, in 
agreem ent with the recent findings of ICescutti fc Chiappinil 
(l2014li . We find that a fraction of 0.1% of all CCSNe which 
end up in this explosion channel provides the best fit. This 
would mean that not all but only a fraction of magnetar 
events which produce the highest magnetic field neutron 
stars are able to eject a main r-process composition of the 
heaviest elements, as discussed above. 

Our conclusion is that both sites acting in combination 
provide the best scenario for understanding [Eu/Fe] obser¬ 
vations throughout galactic history, with typical probabili¬ 
ties for NSM formations and (merging) delay times as well 
as probabilities for Jet-SNe. As a side effect we realized 
that present supernova nucleosynthesis yield predictions, 
based on induced explosions with a single explosion energy 
throughout the whole mass range of progenitor stars, bear 
a number of uncertainties. While apparently too large scat¬ 
ters of alpha/Fe ratios can be obtained in inhomogeneous 
chemical evolution models, when utilizing existing nucle¬ 
osynthesis predictions from artificial explosions with ener¬ 
gies of 1 Bethe, this might not be due to the chemical evo¬ 
lution model. Such deficiencies can be cured by assuming 
larger mixing masses with the ISM for supern ovae explo¬ 
sions llvan de Voort et al.ll201^ . IShen et al.ll2015l l. or the in¬ 
troduction of an artificial floor of abundances based on IMF- 
Integrated yields of CCSNe for metallicities at [Fe/H]= —4, 
but it could in fact just be due the non-existence of self- 
consistent CCSN explosion models. We have shown that an 
explosion energy dependence on the compactness of the Fe- 
core, related to the main sequence mass, could solve this 
problem as well by modifying the nucleosynthesis results. 
Therefore, self-consistent core collapse calculations with ex¬ 
plosion energies varying with progenitor mass and possibly 
other properties like rotation are highly needed. Although 
we have obtained a good accordance with the observed Eu 
abundances, the true origin of r-process elements might still 
require additional insights. The present investigation may 
be used to put constraints on the yields, as well as essential 
properties and occurrence frequencies of sites. 

There exist a number of open questions not addressed in 
the present paper, related (a/b) to production sites and (c) 
to the true chemical evolution of the galaxy. 


(a) As discussed in detail in section 12.2.21 we did not 
include ’’regular” CCSNe from massive stars as contribu¬ 
tors to the main or strong r-process, producing the heav¬ 
iest elements in the Universe. However, as already men- 
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tioned in section 12 .2.21 there exists the chance for a weak 
r-process, producing even Eu in a Honda-style pattern in 
such events. This could provide the correct lower bound of 
[Eu/Fe] in Fig. [5] and would be cons i stent with the recent 
findings of lTsuiimoto fc Shigevam^ (|2014| 'I. 

(b) We did not include the effect of NS-BH merg¬ 
ers in the present paper. They would result in simila r 
ejecta as NS-NS mergers per event JKorobkin et al.ll2012h . 
but their occurrance frequen ci es bear high uncertainties 
( Postoov fc YuneelsorJ 120141 '). iMennekens fc Vanbever^ 
( 20141 ') provide a detailed account of their possible con¬ 
tribution and also discuss their contribution to global 
r-process nucleosynthesis. One major difference with re¬ 
spect to our treatment of NSM in chemical evolution re¬ 
lates to the fact, that (if the black hole formation is not 
causing a hypernova event but rather occurring without 
nucleosynthesis ejecta) only one CCSN is polluting the 
ISM with Fe before the merger event, in comparison to two 
CCSNe. This would lead to a smaller [Fe/H] ratio in the 
ISM which experiences the r-process injection, and just to 
an ’’earlier” appearance of high [Eu/Fe] ratios in galactic 
evolution. If we assume that BH formations are as fre¬ 
quent as supernova explosions, an upper limit of the effect 
would be that all NS-NS mergers are replaced by BH-NS 
mergers, moving the [Eu/Fe] features to lower metallicties 
by a factor of 2. However, the lower main sequence mass 
limit for BH formation is probably of the order 20 Mq, 
and only a small fraction of core collapses end in black 
holes. Therefore, we do not expect that the inclusion of 
NS-BH mergers shifts the entries by more than 0.15 in 
Fig-El This by itself would not be a solution in terms of 
making only compact (i.e. NS-NS and NS-BH) mergers 
responsible for the r-process at very low metallicities. 

(c) There have been suggestions that the Milky Way in 
its present form resulted from merging subsystems with 
a different distribution of masses. Such “dwarf galax¬ 
ies” will experience different star formation rates. It is 
known that different star formation rates can shift the 
relation [X/Fe] as a function of [Fe/H[. If the merg¬ 
ing of such subsystems will be completed at the time 
when type la supernovae start to be important, the re¬ 
lation [X/Fe] =f([Fe/H]) will be uniform at and beyond 
[Fe/H]> —1, but it can be blurred for low metallicities 
between the different systems, possibly leading also to a 
spread of the onset of high [Eu/Fe] ratios a t low metal¬ 
licities (llshimaru. Wanaio fc PrantzosI l2015l l. The result 
depends on the treatment of outflows, should in principle 
be tested in inhomogeneous models, and also already be 
pres ent in the s imula tions of Ivan de Voort et al.l (l2015l l 
and IShen et al.l (l2015li . But it surely requires further in¬ 
vestigations to test fully the impact of NSM on the r- 
process production in the early galaxy. 

Future studies will probably require a distribution of 
delay times for NSM events, a test of the possible contribu¬ 
tions by BH-NS mergers, a better understanding of yields, 
and improvements in understanding mixing processes after 
supernovae explosions and during galactic evolution. Testing 
the full set of element abundances from SNe la and CCSNe 
as well as the two sources discussed above, in combination 
with extended observational data, will provide further clues 
to understanding the evolution of galaxies. 
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